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equation in cylindrical domain by partial boundary data 

O. Yu. Imanuvilov* and M. Yamamoto^ 

O 

(N 
> 

Q I Abstract 

. Let C be a bounded domain with dCl £ C°° and L be a positive number. For a 

\ three dimensional cyHndrical domain Q = Q x {0, L) , we obtain some uniqueness result 

of determining a complex-valued potential for the Schrodinger equation from partial 
. Cauchy data when Dirichlet data vanish on a subboundary (90 \ T) x [0, L] and the 

corresponding Neumann data are observed on F x [0,L], where T is an arbitrary fixed 
■ open set of dCl. 

This article is concerned with the inverse boundary value problem of determination of a 
complex-valued potential for the Schrodinger equation in a cylindrical domain from partial 
boundary data. More precisely the problem is as follow. Let Q = Q x (0, L), where C 
is a bounded domain with dfl G C°°. Let F be an arbitrary, fixed subdomain of dfl. Denote 
Fo = Int{dn \ f ), S = f X [0, L] and Sq = Fq x [0, L]. 

In Q, we consider the Schrodinger equation with some complex-valued potential q: 



>: 



Lg{x,D)u = {A + q)u = in Q. 
Consider the following Dirichlet-to-Neumann map Ag 



Ou 

X : ^I'^of = ^|aQ\So, where Lq{x,D)u = in Q, m|so = 0, u\aQ\^o = f (2) 

with domain 

/^(A,,Eo) = {/e^^W)|supp/cag\So, {f,g)LHdQ) = WgeAf} 

and 

du 

^f = {9\Lg{x,D)u = in Q, u\aQ = 0, —\QQ=g}. 

The problem ([T]) and ([2]) is the generalization of the inverse boundary value problem of 
recovery of the conductivity, which is also known as Calderon's problem (see [2j). It is 
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related to many practical applications, for example, detecting oil or minerals by applying 
voltage and measuring the fluxes on earth's surface. See also Cheney, Issacson and Newell 
[3j for applications to medical imaging of EIT. 

In case when Q is a general domain in with n > 2, Sq = (i.e., the case of full 
Dirichlet-to-Neumann map), the unique recovery of the conductivity was established in [16] 
and [17J in two and three dimensional cases respectively. From the practical point of view, 
the assumption So = , means that one has to set up voltages and measure the fluxes on 
the whole boundary is very restrictive. In practice, this assumption does not often hold, for 
example because the domain Q is extremely large or we can not access to some part of dQ, 
e.g., the domain has cavities located inside. For the inverse boundary value problem with 
such partial Dirichlet-to-Neumann map, we refer to the following works. In [1] Bukhgeim 
and Uhlmann show that if voltages are applied on the boundary OQ^ and the corresponding 
fluxes are measured on some part dQj^ which is approximately equal to dQ \ dQ^, then the 
potential can be uniquely determined. This result and a recent improved result [lOj still 
require the access to the whole boundary dQ. In P], Isakov solves the case where voltage 
applied and current measured on the same set dQ_ provided that subboundary dQ \ dQ^ is 
a part of some sphere or some plane. All the above mentioned papers treat the case where 
the spatial dimension more or equal 3. 

The purpose of this article is to establish the uniqueness with weak constraints on such 
subboundary in the case of three dimensional cylindrical domain Q. Our proof is based on the 
Radon transform and gives the uniqueness in higher dimensional domains with generalized 
geometrical conflgurations (not necessarily cylindrical domains). 

As for related work in slabs, see Ikehata Krupchyk, Lassas and Uhlmann [llj, 
Uhlmann et al [15] . See Novikov [H] for conditional stability results for Calderon's problem. 

For general two dimensional domain, [B] proved the unique recovery of a potential for 
the Schrodinger equation in the case when voltage applied and flux measured both on an 
arbitrary open set of dQ. Thus [6j established the best possible uniqueness in the two 
dimensional case with data A^^s^ deflned by ([2|). Also see [7] which deals with the same 
inverse problem for more general second-order elliptic equations in the two dimensional case 
and |8j improves the result of [B] in terms of regularity assumption of potential for the 
Schrodinger equation. 

We introduce the subset O of domain 

= n\ C/i(f o), Ch{fo) = {x\x = Xx^ + (1 - X)x\ x\x'^ eT^,Xe (0, 1)}. 
We have 

Theorem 1 Letqi,q2 be Lipschitz functions. //Ag^ So = A^j.So o,T^d D{Kq-^ Y.o) C /^(Agj^So); 
then qi = q2 in O X [0, L] . 

From theorem [1] we obtain immediately 

Corollary 2 Let Q be concave nearTo and potentials qi,q2 be Lipschitz functions such that 
Agi,Eo = ^q2,^o ^(Agi,Eo) ^ D{Ag^^Y;o)- Thcn qi = q2 tnQ. 

First we formulate the following Carleman estimate with the linear weight function = xs 
for the Schrodinger operator ([T|). Denote || ■ = || • ||/i'i(Q) + |t||| • ||l2(q). In P the 

following theorem is proved: 
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Theorem 3 Let q G L°°{Q). There exist tq and constant C independent of t such that for 

all T > Tq 

fill 

lke-^||H.^{Q)<C(||(L,(a;,D)w)e-^|U.(Q)+v^||(— e-^)(-,L)|U.(f,)) ^ueHHQ). (3) 
Next we formulate some known results on the generalized Radon transform: 

(7^^/)(u;,p) = I f{x)e^<^^'^>ds, {oo,p) eS'xR, 00^ = {002, -coi). 

The following theorem is proved in [12] 

Theorem 4 Let f be a Lip schitz function with compact support. If {Tlnf){uj,p) = for all 
p> r then supp f C {x G M^||x| < r}. 

Similar to Corollary 2.8 of [4j p. 14 we prove 

Corollary 5 Let f be a Lipschitz function in with compact support and E be a bounded 
convex set in M^. // (7l^f){u,p) = for all lines < u,x >= p which do not intersect E then 

f{x) = yx^E. 

Proof of Theorem [1]. Let point x = {xi,X2,X3) G O x (0,L). Since ChlTo) is the convex 
closed set the point {xi,X2) can be separated from it by some line £. Then the line which 
is parallel to i and pass through {xi,X2) does not intersect Ch{TQ). After possible rotation 
and translation we may assume that Xi = 0, £2 > and axis X2 does not intersect Ch{To). 
Therefore it suffices to prove that 

qi = q2 on Qn {x\x3 G [0,L],Xi = 0}. 

Without loss of generality we may assume that fl C M^. Let m be a smooth function defined 
on such that \m'\ < 1. Denote V = {dx^,dx2)- Consider the eiconal equation 

|V'^f| = 1 inQ, ^f|x2=o = m. (4) 

This equation can be integrated by the method of characteristics. The solutions, as long as 
they exist, have the following form 

\l/(xo + ta{xo)ei + t/3(xo)e2) = m{xo) + t Vxq G {x'\x2 = 0}, t > 0, (5) 



where a{xo) = m'{xo), /3{xo) =a/1 — a^(xo). 

Next construct the function \E' more explicitly using the implicit function theorem. Con- 
sider the following mapping: 

F{y) = y2a{yi)ei + y2(3{yi)e2 + (^^^^ , y = {yi, 2/2)- 

Assume that 

a{0) = m'(0) = 0. (6) 



3 



Then 

F(0,t) = (0,t) (7) 

and 



In particular 



F'{0,y2)-- 

As long as the function 1 + y2a'{0) is positive, there exists the inverse matrix 



l + l/2«'(0) 
1 







(F')-^(0,2/2)= J. (8) 

Let i^' be a positive number such that 

n n {x'\xi = 0} C {x'\xi = 0,0 < X2 < K}. 

By ([8]) there exists e{K) for any a e X = {(f) e C^{-1, 1)|0(O) = 0} such that ||a||c2[_i,i] < 
e{K) there exists 6 > such that on the set [—5,5] x [0,2i^] the matrix {F')^^ is correctly 
defined 



Then by ([7j) and the implicit function theorem there exist 5 > such that the mapping 
x' — 7- y{x') is correctly defined on 2) = [—5, 5] x [0, i^T] and the derivative of this mapping 
given by formula: 

^ = {Fr\y{x')). (9) 

Differentiating the first columns on both sides of the matrix equation ^ with respect to Xi 
and using ([6]) we have 



(0' ^2) = -JT-T-TTTT^^ n + /1 , ^./^n^^2 • (10) 



; ~ "(i+i/2«'(o))3 vo; ^ (i+y2«'(o))2 v-i/2(«'(o))^ 

Then the function \E' can be determined by formula 

^(x') =mo(yi(x')) +?/2(a;')- 
The short computations and ([HD, ffTU]l imply 

Av^(0,X2) =m[;(0)(9.,yi)2(0,X2)+m[,(0)(a^^l/i)(0,X2) + 9,^1/2(0, X2) - °^ ^ 



;i + Z/2a'(0))^ 



y2{a'{0)f _ a'{0) 1/2^(0))' _ a^(0) 

+ (l + 2/2a'(0))2 " (l + y2«'(0))2 + (l + y2a'(0))2 " (T+^^^W ^^^^ 
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Let ao{x') be a function such that 

2(V'^, V ao) + A^./^ao = in D, (12) 
and a{x') be a smooth function such that 

(V^, V'a) = in O. (13) 
Next we construct the functions ao and a. 

In order to construct the function a{x') we take a smooth function r 

reC^{-e,e), (14) 

where e is a small parameter and set: 



a(x') = r(xo) on the hne {x' E M.'^\x' = (^q^ + ta{xo)ei + 



t/3{xo)e2,t>0}. (15) 



We claim that for the function a defined by these formula we have f[T3|) . Set Vi = a(xo)ei + 
f3{xo)e2,V2 = a{xo)ei — f3{xo)e2- Then by dSD If^i/i^l = 1- Since |V'\1/| = 1 we have that the 
vector Vi is parallel to V'\l'. Hence vectors Vj are orthogonal. So 



Therefore 



d^^'^ = 0. 



d^,^a = |V'^|5v'* a = \V-^\d^,a = 0. 

|V'*| 



Hence the formula (fT3|) is proved. 

We integrate equation (fT2l) by the characteristic method. In particular using ( [TT]) we 
have 

ao(0,X2) = e"2^° (i+«2-'(0))'^?'2 = -^«n(l+X2a'(0)) _ (16) 

Next we construct the complex geometric optics solution ui{x,t) for the Schrodinger 
operator with the potential qi. For the principal term of complex geometric optics solution 
we set 

U = e(^+^)(^«+^*(^'»aao. (17) 

The set O is closed and the axis X2 does not intersect this set. So there exists a neighborhood 
of the set {x'|x2 G [0, K]} such that it does not intersect O. Thanks to (IT^ and (ITSl) . choosing 
a positive parameter e sufficiently small, we obtain 

[/|so = 0. (18) 

The simple computations imply 

L,, (x, D)U = (r + N)\V{x3 + i^), V(x3 + 2^))f/ 
+ (r + Ar)(2(V'^, V'ao) + A,,*ao)ae("+^)("^+'*("'» 
+e(r+iv)(-3+^*(x.'))^^,(^^^) + q^U = e(^+^)(^'^+**(^''»A,.(aoa) + q,U. (19) 
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Observe that the functions (e*-^"*"^-'^^^'^***^^'-'^ A^/(aoa) +gif/)e~™^ are uniformly bounded 
in T in norm of the space L'^{Q). Consequently using the results of [1] we construct the last 
term in complex geometric optics solution- the function Ucor{-, t) such that 

Lg^{x,D){e'^'^Ucor) = -Lg^{x,D)U in Q, Ucorlsi = (20) 

and ^ 

WucotWl^q) = 0{-) asr^+cx). (21) 

Hence, by ( 1T8|) . ( I20l) and ( 12T]) we have the representation 

u^ = U + e^^'OL^Q){^) as r ^ +cx), MiIso = 0. (22) 

(By Ol^(^q) we mean any function f{-,T) such that ||/(-, t)||l2(q) = 0(1) as r — > +oo. ) 
Similarly we set 

V = e-"("^+'*("'))ao, y|so = 0. (23) 

We multiply any smooth function qq satisfying (fT2!) with a solution of equation (fT3|) which 
is supported around the ray {x|x2 > 0,xi = 0} and is equal to 1 on this ray. Hence we can 
assume that the support of function gq is concentrated around this ray and (fT6|) holds true. 

Since the Dirichlet-to-Neumann maps of the Schrodinger equations with potentials qi, q2 
are the same there exists a solution to the following boundary value problem 

Lqoix, D)u2 = in Q, ui = U2 on dQ, and ^— = — — on dQ \ Sq. (24) 

Oh' ou 

Setting u = Ui — U2 and using f l2^ we have 

du 

Lg^{x, D)u = (gi - q2)ui in Q, u\dQ = ^|9q\So = 0. (25) 

Applying to equation ( !25|) the Carleman estimate ([3]) we have that there exist constants 
C and To independent of r such that 

lke-^1li/i-(Q) < C Vr>ro. (26) 

Then taking the scalar product in L'^{Q) of equation ( |2^ with V, and using ( l26ll . ( |23l) . 
fl22|) we have 

I (gi — q2)uiVdx = f uLg.^{x, D)Vdx = 0(— ) as r — )■ +oo. 
Jq Jn 

This equality and ( I2T1) imply 

/ (gi - g2)e'^("3+**("'»aagrfx = O(-) as r ^ +oo. 
Jq ^ 

Therefore 

{qi - g2)e'^("^+'*("'))aa2rfx = 0. (27) 

Q 
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Setting pn{x') = J^iqi — q2)e^^^dx^ and using ( !T6|) we obtain from (j2 



K 



PNe'^''^dx2 = 0. (28) 

Indeed, let r = be a function such that it is equal l/2h on the segment [—h,h] and 
zero otherwise. Denote the solution to equation (fT3|) given by (fT5|) with the initial condition 
r/j as a{h). By (ITSl) the function ao(/i) given by formula 

x'^n^, 



a{h) 



where II/j = {x' G M^|x2 G [0, -fC] — /i + ^p^X2 < Xi < /i + ^|j^X2}. Therefore for any fixed 
X2 from the segment [0,fr] the function a{h) equals on the segment [—h + ^737^X2,/?- + 



^ on the bv^g,^^^-^^- L '"I /3(-h) 

m 



^^X2]- By ([6]) the length of this segment is 2h + 2a'{0)x2h + o{h). We rewrite f l27|) as 



= / (gi-g2)e^(^^+^*(^'»a(/i)a^rfx = ^ / p^e^^^^^'^a^^x' = / / "° rfxiC^X2. 



2^ /n, 7o 2/1 



Applying (fT6l) and using the assumption on regularity of potentials qj we have 

rK rh+a'io)Hx2 ^2^m^ fK i-h+a'mhx, p^(0, ^,)e^A.^2 + o(l) , 
= / / dxidx2+o{l) = / — — — — - — r — dxidx2+o{l). 



Jo J-h-a'{0)hx2 2/1 Jq J-h~a'{0)hx2 2/l(l + «'(0)X2) 

This proves (|28|). 
Next we claim 

P7ve"*^^'rfx2 = 0. (29) 

The proof of ( l29l) is exactly the same as the proof of equality ( l28i) . The only difference 
is that instead of the function \1/ one has to use the function defined by 

IV'^I = 1 in n, ^1^2=0 = rn, (30) 

\E'(xo — ta;(xo)ei + t/3(xo)e2) = m{xo) — t Wxq G {x'|x2 = 0} and Vt > 0, (31) 

where a{xo) = m'(xo), /3(xo) =a/1 — a^(xo). 

From fl28l) and fl29|) setting = —i'-f where 7 is the real parameter we have 

n^{p-^J){u,p) = y{u,p) e S'^ xR^ such that {x\ < u,x >= p} nch{n\o) = ^. 

Applying corollary [5] we have that 

p-i^{x') = W eO and V7 G 

Therefore for any fixed x' E O and any 7 

r-T 

(qi - q2)ix',X3)e'^''''dx3 = 0. 



This equality implies immediately that the function X3 — )■ (gi — q'2)(x',X3) on the segment 
[0, L] is orthogonal to any polynomials. Therefore (gi — q2){x)\ox[o,L] = 0. The proof of the 
theorem is complete. 
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